How universal is the one-particle Green's function of a Luttinger liquid? 
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The one-particle Green's function of the Tomonaga-Luttinger model for one-dimensional interacting 
Fermions is discussed. Far away from the origin of the plane of space-time coordinates the function 
falls off like a power law. The exponent depends on the direction within the plane. For a certain 
form of the interaction potential or within an approximated cut-off procedure the different exponents 
only depend on the strength of the interaction at zero momentum and can be expressed in terms of 
the Luttinger liquid parameters K p and K a of the model at hand. For a more general interaction 
and directions which are determined by the charge velocity v p and spin velocity v a the exponents 
also depend on the "smoothness" of the interaction at zero momentum and the asymptotic behavior 
of the Green's function is not given by the Luttinger liquid parameters alone. This shows that the 
physics of large space-time distances in Luttinger liquids is less universal than is widely believed. 

PACS numbers: 71.10.-w, 71.10.Pm 
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Over the past 20 years it has been shown that the low 
energy physics of a variety of models of one-dimensional 
(ID) correlated electrons can be described by the Lut- 
tinger liquid (LL) phenomenology Because of the 
progress in the experimental realization of quasi ID sys- 
tems and speculations about possible LL behavior in the 
normal state of the high-temperature superconductors 
LL phenomenology has lately attracted considerable at- 
tention. In Fermi liquids the low energy physics and low 
temperature thermodynamics is dominated by the exci- 
tation of quasi-particles which are in a one-to-one cor- 
respondence to the particle-hole excitations of the non- 
interacting system. The elementary excitations in LL's 
are of collective bosonic nature. In LL's the discontinu- 
ity of the momentum distribution function n(k) at the 
Fermi wave vector fej? vanishes and the density of states 
of the one-particle Green's function at the Fermi energy 
is zero. Close to the Fermi surface both functions are 
dominated by power law behavior with exponents which 
are given by the so called LL parameters K p and K a . 
According to LL phenomenology thermodynamic quan- 
tities (e. g. the compressibility and spin susceptibility) 
and the non-analytic behavior of correlation functions 
can be expressed in terms of these two parameters and 
the velocities v p and v a of charge and spin excitations. 
By calculating the one-particle Green's function of the 
Tomonaga-Luttinger (TL) model ||Q we will show that 
the asymptotic behavior of the Green's function of a gen- 
eral LL at large space-time distances is nonetheless less 
universal than is widely believed. The exponents of the 
algebraic decay in the directions within the x-Vpt plane 
which are determined by the velocities v p and v a are not 
given by the LL parameters of the model alone. Here vf 
denotes the Fermi velocity. We will furthermore discuss 
the implications for the spectral function. 

The TL model is a continuum model of interacting 
ID electrons. The linearization of the electron disper- 
sion around the two Fermi points and neglecting the 
backscattering processes between electrons makes it fea- 
sible to determine the spectrum of the Hamiltonian and 



calculate all correlation functions. Following Luttinger 
Q we introduce right- (a = +) and left-moving (a = — ) 
Fermions with spin s, creation operators at , disper- 



sion £ a (k) = 

PaAl) = v < " 



avF(k — akp), density operators (q =/= 0) 
k a k,a.s a k+ q ,a,s' an d particle number opera- 
tors Uk,a,s = a \ a s a k as- To simplify the mathematical 
treatment Luttinger added an infinite filled Fermi sea to 
the ground state pj. The Hamiltonian for a system of 
length L is given by 



q^0 



[52,11 (e)*^' + 92,±(q)S s -s'] p+, s (q)pl_ s ,(q). (1) 



Here we use "g-ology" notation M and (. . .) denotes the 
(non-interacting) ground state expectation value. Con- 
trary to many authors we keep the explicit q dependence 
of the coupling functions. We assume that the Fourier 
transforms gi, K {q) (i = 2,4; n =||,_L) of the two-particle 
interaction have only contributions for q^q c <SC fc^ with 
an interaction cut-off q c . At no stage of the discussion 
it will be necessary to introduce any further cut-offs "by 
hand" despite the infinite (filled) Fermi sea at negative 
energies. The model only belongs to the LL universality 
class if g2, K {q — 0) is finite for k —\\ and _L and at least 
one of the two coupling constants is non-zero. Thus we 
restrict ourselves to these kind of interactions. 

Bosonization of the Hamiltonian and a canonical trans- 
formation leads to || 



H 



2J 



q^0 v—p,a 



(2) 



with bosonic operators (3 v {qy describing charge (y = p) 
and spin (y — a) excitations (spin-charge separation). 
The energies e q , v are given by 



1 



= V F 




g±Aq) 

irvp 



TTVp 



= v„{q), (3) 



where we have introduced the renormalized charge and 
spin density velocities v v (q) and interactions gi :P / a (q) = 
[^^(g) ± gi,±(q)] /2. The one-particle Green's function 
iG< s (x,t) = (-0q, s (O, 0)7/) a s (a;, t)), which after a dou- 
ble Fourier transformation leads to the spectral func- 
tion /?< s (fc,o>) relevant for photoemission experiments, 
can be calculated using the bosonization of the fermion 

(TJ. In the ther- 



fields ipi !S (x) = ^E fc e 



-ikx „1 



k,a,s 



modynamic limit and at zero temperature we obtain 
G<(x,t) = [G<]°{x,t)exp{F{x,t)} with 



-iq(>-u„(q)t) 



u=p,a J ° y 

+ 2-y u (q) [cos^e^-W*-!]} 
and the non-interacting Green's function 



-iq(x—Vj?t) 



[G<?(x,t) 



1 



2tt x — vpt — iO 



(4) 



(5) 



Because the Green's function is the same for both spin 
directions we have suppressed the spin index s. By leav- 
ing out irrelevant particle number contributions in the 
Hamiltonian we effectively shifted the energy scale so 
that the chemical potential is zero, jviq) in Eq. (^) is 
given by the eigenvectors of the canonical transformation 
and can be written as jv(q) = [K v (q) + 1/K v (q) — 2] /4, 
with 



K v {q) 



g^u{q)/{nv F ) - g2,v{q)/(irvF) 
(q)/ (ttv f ) + g2, v {q)l (ttv f ) ' 




(6) 



Due to the interaction cut-off the momentum integral in 
Eq. (Q) is regular in the ultraviolet limit. The q — > 
limits of the four functions v v (q) and K v (q) define the 
two velocities and two LL parameters (v p ,v a , K p , K a ). 

Roughly speaking the non-analyticities in spectral 
functions at low energies and small momenta are given 
by the behavior of G^(x,i) at large arguments. To 
determine the non-analytic behavior of the momen- 
tum distribution n + (k) = J^° oo dxe~ ikx iG+(x,0) and 
the momentum integrated spectral function p^(uj) = 
dte tuJt iG< (0, t) I (27r) it is sufficient to discuss the be- 
havior of G+ (x, t) along the x and VFt axis. To emphasize 
the difference to the general case of non-vanishing x and 
v F t we first focus on these two cases. With a little less 
mathematical rigor they have already been discussed sev- 
eral times ||. Using integration by parts Q in the first 
derivative of F(x,0) and F(0,t) we obtain the leading 
behavior F(x,0) ~ — aln|x| and F(Q,t) ~ —a\n\vpt\, 
with a — 7 P (0) + 7,7(0). In G+(x, t) this leads to power 
law behavior along the x and v Ft axis with an exponent 
which only depends on the strength of the interaction 



at vanishing momentum and thus the LL parameter K p 
and K a of the model at hand. This is in accordance with 
the LL phenomenology. From general theorems about 
Fourier transforms the by now well known LL behav- 
ior of the momentum distribution function follows for 
\k-k F \/q c ^0 



\k — kp\ a sign(k — k F ) 
.(k)~{ {k-k F )\n\k-k F \ 
(k - k F ) 



for < a < 1 
for a = 1 
for a > 1. 

(7) 



For a > 1 the leading behavior of n + {k) is dominated by 
a linear term but a higher derivative still diverges at k F - 
Due to the different analytic properties of F(0,t) (ana- 
lytic in the upper half of the complex t plane) and F(x, 0) 
(not analytic in either the upper or the lower half of the 
complex x plane) P+(oj) is dominated by a non-analytic 
power law behavior even for a > 1 Q. Thus anomalous 
dimensions a > 1 should in principle be observable in 
momentum integrated photoemission spectra J^]. In the 
literature the mathematical reason for this important dif- 
ference in the behavior of n+(k) and p+(w) has thus far 
not been properly pointed out. For oj/(v F q c ) — > we 
obtain 



@{-u>){-u>r 

e(-w)(-w) Q ln|w| 



for a # IN 

for a € IN. 1 ' 



The prefactors and the range over which the leading be- 
havior given in Eqs. (0) and (|§|) can be observed depend 
on the interaction at all q. As a consequence numeri- 
cally calculated curves for different interaction potentials 
gi.n{q) but the same anomalous dimension might appear 
quite different 

To find a more explicit form of the Green's function 
the q integral in Eq. (|J) has often been evaluated after 
replacing v v (q) — ► i>„(0), lv(q) — > 7^(0) and multiplying 
the integrand by a factor exp (— qA) |||ll) . The q integral 
can then be performed and one obtains 



[G<} A (x,t) 



— e 



7(2*0 



x — Vpt 

A 2 



iO 



n 

y= P: a 



X 



Vpt — iA 



(x — v v t — iA) (x + v v t + iA) 



x — v v t 

7x/ 2 



iA 



1/2 



(9) 



For the Hamiltonian Eq. ([!]) there exists no special in- 
teraction potential so that this approximation becomes 
exact. As in the discussion of the general case below 
we will transform onto new variables s = x — ct and 
s' = x + ct with an arbitrary velocity c and discuss the 
behavior of [G+J^ (x[s, s'],t[s, s']) for large s with a fixed 
s' and vice versa and for different values of c. For all 
velocities c but v p and v a , [G^]^ (x[s, s'}, t[s, s']) falls off 
like s - ( 1+Q ) and s'~( 1+a \ This is the behavior we al- 
ready found for vpt = or x = 0. For c = v p the Green's 
function falls off like 



2 



-(l+7<r+7p/2) 



[G<] A (x[ s , s '],t[ s , s ']) , 

[G<U(x[ S , S '],i[ S , S ']) ^ S '-(l/2+7.+7 P /2). 



(10) 
(11) 



The behavior of the Green's function for c = v a fol- 
lows from Eqs. (|Io| ) and ( |ll| ) by interchanging 7 p and 
7 CT . Within the above approximation the exponents of 
the asymptotic behavior for all c and thus for all direc- 
tions within the x-vpt plane can be expressed in terms 
of K p and K a . The resulting spectral function dis- 
plays power law singularities (bounded non-analyticitics 
in case the 7^ are too large) at u) — ±t>„(fc — k F ) for all 
A(-HM>0§. 

Without using any approximations it has been shown 
that Eq. (^) gives the correct asymptotic behavior of 
Gl(x,t) for a box potential g^ K (q) = g^Sfe - \q\) 
|[10|| . In this case the momentum range over which 
p+(k,uj) is dominated by two power law singularities 
at the charge and spin excitation energies is limited to 
< (—k + kp)/q c < 1. For other momenta further non- 
analyticities occur |Tc|| . 

Next we will show that [G^]a(x, t) does not give the 
correct asymptotic behavior for an arbitrary shape of the 
interaction. After transforming onto s and s' the Green's 
function is given by 

G< (x[s, s'],t[s, s'}) = [G<}° (x[s, s'],t[s, s'}) e^ s > s '\ (12) 

where F(s.s') follows from Eq. (0) and (x, i) — > {s,s'). 
As an example we will discuss the behavior of F(s, s') for 
large s' and fixed s in more detail and only present the 
results for the other case. We first take the derivative 
with respect to s'. This gives 



dF(s,s') 1 



xe 



u—p.a 

3 - 

e 



dq 



4< 



v v {q) 



ds 1 2 

- lg l(i+-M) Se - jg l(i--M) s ' 



1 



c / 

.1 L + MoT } ^ (i-a^i). e i,j (i+^m) 



(13) 



To simplify the discussion we assume that v v (q) are 
monotonic functions. Using integration by parts and the 
method of stationary phase in the asymptotic expansion 
of the integral we find for all velocities c but v p (0), v a (0), 
and vp that F(s,s') goes like — aln(s') and thus 



at q = 0. At the stationary point also the prefactor 
[1 — Vp(q)/v p (0)] vanishes and thus we have to general- 
ize the method of stationary phase. The details of this 
generalization will be given elsewhere and here we will 
only present the results. The leading contribution of 
the fourth term of Eq. (|l3|) to the large s' behavior is 
— ^j- - 1 +1 jj, where p p is the smallest number n € INU{oo} 

with [u P ]^"' ) (0) ^ 0, where [v P ]^ (q) denotes the n-th 
derivative. According to the definition of v p (q) Eq. (||) 
Pp is a measure of the "smoothness" of the interaction at 
vanishing momentum. Evaluating the first term in Eq. 
( |l3| ) in the same way and using integration by parts in 
the other terms leads to 



G< (s,s')~s 



/- (1/2+7, +7 P /2+l/ [2p p +2] + 7p / [2p p +2] ) 



(15) 



This result is different from Eq. (11) and more impor- 
tantly it shows that the asymptotic behavior cannot be 
obtained from the LL parameters K p and K a of the TL 
model alone. This is in contrast to the wide spread belief 
that Eq. (|^) displays the asymptotic behavior of all mod- 
els which belong to the LL universality class provided the 
exponents are expressed in terms of the LL parameters of 
the specific model. The behavior Eq. ([ll]) is only recov- 
ered if p p — oo, i. e. if all derivatives of v p (q) and thus all 
derivatives of the interaction potential vanish at q = 0. 
For this reason Eq. (^|) gives the correct leading behavior 
in case of a box potential. For the large s behavior and 
c = v p (0) we obtain 



G< 



l+7„+7p/2+7 P /[2p P +2]) 



(16) 



For c = v a (0) we have to analyze the Green's function 
following the same route and obtain 

G< (S, S') ~ s -(l+7 P +7./2+7./[2p„+2]) ; 

G<(s,s') ~ s '-(l/2+7p+7^/2+l/[2p„+2]+ 7 „/[2p^+2])^ ^ 

where p a is defined in analogy to p p . 

The case c = vp has so far been excluded. If vp ^ 
v p (0) and vp ^ v a (0) the behavior of the Green's func- 
tion is given by Eq. fl4|). In the relevant case of a spin 
independent interaction with v a (q) = vp and 7 CT = 
the Green's function is given by Eqs. ( p^?| ) and (18) with 
Pa = oo. 

To gain more insight into the asymptotic behavior of 
the Green's function we have evaluated the momentum 
integral in F(s, s') numerically. As an illustration of the 
new predictions Eqs. (|15|)-(|l8|) we present Re{F(s,s')} 



for s = as a function of s' and c = v p (0) on a log- linear 
scale in Fig. |. The curves have been calculated for a 
spin independent interaction with 34,11(9) = gi,±(q) = 
52,||(9) = 92.±{q) = gexp{(-q/q c ) p "}, 2g/{-Kv F ) = 3, i. 



G< (s, s') = G< (x[s, s'],t[s, s']) - s' _(1+q) . (14) e. v p (0)/v F = 2, 7p = 1/8 and different p p . From Eqs 



For large s and fixed s' we obtain the same behavior. 

If c = v p (Q) the phase q fl — V o (q)/v o (0)] in the 
first and fourth term of Eq. ( jl3| ) becomes stationary 



dl2] ) and (jl5|) we expect 



F(s,s>) 



It 
2 



1 



7p 



2p p + 2 2 Pp + 2 



In(s') (19) 



3 



Fits of the data for s' > 10 3 reproduce the expected 
prefactors of the In (s') within a relative error of less 
than 0.02%. For increasing s we have to go to larger 




FIG. 2. Re {F(s,s')} as a function of s'q c for different sq c . 
For details see the text. 

s' to find the asymptotic behavior Eq. (|l9|). In an in- 
termediate regime of s' , which increases with increas- 
ing s, Re {F(s, s')} displays the behavior of the p p = oo 
case. This is illustrated in Fig. |. It shows Re{F(s,s')} 
as a function of s' on a log-linear scale for the above 
form of coupling functions, c = v p (0) and different s. 
The parameters are 2g/(wvp) = 5, i. e. v p (0)/vp = v6, 
7 P = 0.2144 and p p — 2. A fit for sq c = again repro- 
duces the expected behavior Eq. ( fL9| ) with a very high 
accuracy. The data for sq c = 10 only show this behavior 
for s'q c > 10 5 . For sq c = 1000 a fit for 10 6 < s'q c < 10 7 
gives 0.3934 In (s') which is very close to 0.3929 In (s'), the 
expected behavior for p p = oo. Similar to the sq c = 100 



curve for arguments between 10 6 and 10 7 , the curve for 
sq c = 1000 shows a cross over at very large s'q c and the 
prefactor of the logarithmic term is again given by the 
p p = 2 value 0.1904. 

For the approximated Green's function Eq. (^) and 
the box potential the two-dimensional Fourier transfor- 
mation which leads to the momentum resolved spectral 
function can be performed analytically pUfiOf . Thus far 
we have not succeeded in calculating p< s (k,uj) for an in- 
teraction with p v < oo. From the Fourier transformation 
of the approximated Green's function it is known, that 
the exponents of the algebraic decay along the special 
directions (c = v v ) determine the exponents of the non- 
analyticities at u) = ±tv(fc — kp). As we have shown 
above the exponents of the algebraic decay of Gf (x, t) 
along the special directions are different from the ex- 
ponents of [G+]ji_(x,t), thus we have no reason to be- 
lieve that p^(k,uj) shows power law singularities with 
the same exponents as \pz ]aQz, uS). It is not even obvi- 
ous that the exact spectral function shows power law sin- 
gularities at all. Certainly we expect the two peak struc- 
ture of spin and charge excitations, but it is not clear 
if, for any non- vanishing k — kp, these peaks are given 
by algebraic singularities. The extended region in which 
G+(x,t) displays the asymptotic behavior of [G^]^i(ar, t) 
on the other hand indicates that the resulting spectral 
functions might look very similar at least for very small 
\k — kp\. For \k — kp\ — ► we expect to find grow- 
ing regions in which the exact spectral function resem- 
bles the power law behavior of the approximation (with 
the same exponents as the approximation) up to energies 
very close to ±v v (k — kp) but not exactly at these ener- 
gies. A comparison of broadened spectral functions for a 
finite system, i. e. of spectral functions with no "real" al- 
gebraic singularities, for a box and Gaussian potential is 
presented in Ref. |8| and indeed shows a prominent simi- 
larity between both spectra. A more detailed comparison 
will be given elsewhere. 

The results presented here are important for the in- 
terpretation of numerically calculated Green's functions 
and spectra of microscopic models. They show that it is 
impossible to determine the LL parameters of the consid- 
ered model from the asymptotic behavior of the Green's 
function within the special directions. The numerical 
evaluation of the asymptotic behavior of the Green's 
function gives on the other hand a possibility to confirm 
our predictions. Furthermore it should be possible to 
confirm the predictions by analyzing the finite size scal- 
ing of the weight of single peaks in the spectral function 
of a microscopic model of finite length. For energies close 
to but different from ±v u (k — kp) we expect to find power 
law behavior with exponents given by the LL paramters 
as in the approximation. Not so for the scaling of the 
peaks exactly at ±i>„(A; — kp). Our results have also con- 
sequences for the comparison of angle resolved and angle 
integrated spectral functions which have been measured 
by high resolution photoemission Q . 
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